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Abstract 

Q 2 evolution of structure functions in the nucleon and nuclei is investigated by 
using usual DGLAP equations and parton-recombination equations. Calculated 
results for proton's F2 and for the ratio F^/F® are compared with various 
experimental data. Furthermore, we study nuclear dependence of Q 2 evolution 
in tin and carbon nuclei: dfF^/F^/dfln Q 2 ] 7^ 0, which was found by NMC. 

1. Introduction 

Internal structure of the nucleon can be investigated in high-energy lepton- 
nucleon scattering. Measured structure functions depend on two variables, Q 2 = 
—q 2 and x = Q 2 /2P ■ q, where q is the four-momentum transfer and P is the 
nucleon momentum. Their Q 2 dependence can be calculated within perturbative 
QCD. An intuitive way of describing the Q 2 variation is to use the DGLAP 
equations Q. They are often used in experimental analysis and also theoretical 
calculations, so it is important to investigate the numerical solution of these 
equations. 

In addition, Q 2 dependence of nuclear structure functions becomes increas- 
ingly interesting. It is because the NMC (New Muon Collaboration) showed Q 2 
variations of the ratio F 2 A /F 2 D with reasonably good accuracy ||. Furthermore, 
it is found recently that there exist significant differences between tin and carbon 
Q 2 variations, d[F2 n / F2}/ d\hi Q 2 ] ^ 0. It is the first indication of nuclear effects 
in the Q 2 evolution of F2 and is worth investigating theoretically. However, it is 
not obvious whether DGLAP could be applied to the nuclear case. In particular, 
the longitudinal localization size of a parton with momentum fraction x could 
exceed an average nucleon separation in a nucleus if x is small (x < 0.1). In this 
case, partons in different nucleons could interact and the interactions are called 
parton recombinations (PR). It is considered that their contributions enter into 
the evolution, and modified Q 2 evolution equations are proposed in Ref. ||. 

The purpose of our study is to investigate the numerical solution of these 
Q 2 evolution equations. We calculate the Q 2 variation of proton's F2 and the 
ratio F2 a /F 2 D . Nuclear dependence of Q 2 evolution d[F2 n / F 2 c ]/d[ln Q 2 ] is also 
calculated in order to understand the NMC measurements. 
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2. Numerical solution 

The DGLAP and PR evolution equations are given by following integrodif- 
ferential equations: 
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where the variable t is defined by t = — (2//3q) ln[a s (Q 2 )/a s (Qo)]- The func- 
tions qi(x,t) and g(x,t) are flavor-i quark and gluon distributions, respectively. 
P Pi p.(z) are called splitting functions which determine the probability that a 
parton pj with the nucleon's momentum fraction y splits into a parton pi with 
the momentum fraction x and another parton. In the PR evolution case, there 
are additional higher-twist terms which are proportional to 1/Q 2 . These terms 
are also proportional to A 1 / 3 because recombination effect is proportional to the 
magnitude of parton overlap in the longitudinal direction. Furthermore, there 
is an extra evolution equation for a higher-dimensional gluon distribution Ght- 
Explicit expressions of recombination contributions are found in Ref. [||, £|. 

As a numerical method for solving the evolution equations, we have been 
studying a Laguerre-polynomial method ||. It is very efficient by considering 
computing time and numerical accuracy. However, this method has some dif- 
ficulties, e.g. in handling the non-linear recombination terms. Therefore, we 
decide to employ a brute-force method M . In this method, the variables x and t 
are divided into N x and Nt steps and integration and differentiation are defined 
by df(x)/dx = [f(x m+l ) - f(x m )]/Ax m and fdxf(x) = J2m=i^ x m f{x m ). If 
initial distributions are given, we can solve the Q 2 evolution equations step by 
step. This is the simplest method for solving the integrodifferential equations, 
but the number of steps N x and Nt have to be large enough to get accurate re- 
sults. Furthermore, the small x region becomes increasingly important with the 
development of high-energy accelerators such as HERA. So it is necessary to have 
a good accuracy at small x as small as 10 -5 . In order to satisfy this condition, 
the logarithmic-x step A(log\ox) = \logiox m i n \/N x is taken in our analysis. 
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It is useful to have a computer program for solving the evolution equations 
accurately because they are frequently used in theoretical and experimental stud- 
ies. Therefore, we provide a FORTRAN program BF1 M which is available at 
the following homepage 

htt p : / / www . cc . saga-u . ac.jp/saga-u/ r iko / physics) /quantuml/structure.html, 

File name : bfl.fort77.gz. 

This is a very useful program for studying spin-independent structure functions 
in the nucleon and nuclei. 



3. Results 

Using the program BF1, we calculate 
various Q 2 evolution. First, we calcu- 
late Q 2 evolution of proton's F% struc- 
ture function by using leading-order (LO) 10 - 
DGLAP, next-to-leading-order (NLO) 
DGLAP, and parton-recombination (PR) ° 
evolution equations. As initial distribu- vf 
tions, we choose MRS(G) distributions 
which are given at 4 GeV 2 . The results 
are shown with various experimental data 10 ' 

in Fig. 1. In this figure, our results agree 10-' 10' 10' io z 10 1 10 4 

with the data very well except for the data QZ (Qev2) 

at small x and at small Q 2 where pertur- 

bative QCD would not work. However, we Figure 1: Q 2 dependence of proton's 
can not test the recombinations because F 2 . 
their effects in the nucleon are very small. 

Next, we investigate nuclear cases. In order to calculate Q 2 evolution of 
nuclear F% structure functions, we need to have input parton distributions in 
nuclei at certain Q 2 . At this stage, there are various models || which can explain 
x dependence of measured ratios F^/F^- In such models, we have a hybrid model 
with parton-recombination and Q 2 rescaling mechanisms M, ||] . According to this 
model, nuclear parton distributions are calculated at Qq=0.8 GeV 2 . We evolve 
these initial distributions to those at larger Q 2 

In Fig. 2, we show the results of Fg a /F£ at x=0.0085 with NMC exper- 
imental data [||. The solid, dashed, and dot-dashed curves are obtained by 
LO-DGLAP, NLO-DGLAP, and PR evolution equations respectively with A=0.2 
GeV and Nf=3. As shown by the figure, NLO and recombination contributions 
to the ratio are conspicuous at small x. If we evolve F 2 from Ql=0.8 GeV 2 , the 
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Figure 3: Nuclear dependence of Q 2 
Figure 2: Q 2 evolution of Ff' a /F<p. evolution. 

recombination effect is larger than the NLO one. It is interesting to find such 
large recombination contributions. However, it is obvious that the recombination 
cannot be tested at this stage because we do not have data at large Q 2 with small 
x. In order to investigate the details of recombination, we need data in the wide 
Q 2 region at small x. 

We also investigate the nuclear dependence of Q 2 evolution in tin and car- 
bon nuclei |J. Calculated results for <9[F 2 /i^J/c^lnQ 2 ] are shown at Q 2 =5 
GeV 2 together with preliminary NMC data fl(| in Fig. 3. The dotted, solid, 
and dashed curves correspond to LO-DGLAP, NLO-DGLAP, and PR evolution 
results respectively. The DGLAP evolution curves agree roughly with the ex- 
perimental tendency. It is interesting to find that the PR results disagree with 
experimental data. Because of the significant discrepancy from the data, large 
parton-recombination contributions could be ruled out. However, it does not 
mean that the parton-recombination mechanism itself is in danger. Actually, 
there is an essential parameter Kht, which determines how large the higher- 
dimensional gluon distribution is (xGht{x, Qq) = KnT[xg(x, Q 2 ,)] 2 )- The mag- 
nitude of Kht is unknown at this stage and we choose Kht = 1-68 from Ref. 
H . In order to discuss the validity of the PR evolution, the constant Kht must 
be evaluated theoretically. It is encouraging to study the details of the recom- 
bination mechanism in comparison with the NMC data. On the other hand, 
proposed HERA nuclear program should be able to clarify this issue by taking 
small x (~ 10 -4 ) data in Fig. 3. 

Finally, we comment on a spin-dependent case. Because NLO spin-dependent 
splitting functions are recently evaluated, we investigate Q 2 evolution of spin- 
dependent structure functions with the NLO contributions. This study is still in 
progress, and it is partly discussed in Ref. [jO]] . 
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4. Conclusions 



We investigate numerical solution of DGLAP and PR equations. We provide 
a FORTRAN program BF1 for solving these equations in a brute-force method. 
Using this program, we calculate the Q 2 evolution of proton's F2 and the ratio 
i*2 '/F 2 L '• Our results are consistent with experimental data. However, parton 
recombinations in the nucleon cannot be found because their effects are very 
small. Their contributions cannot be tested at this stage even in the nuclear 
case. In order to investigate the details of the recombination in nuclei, we need 
data in the wide Q 2 region at small x. We also find the nuclear dependence of Q 2 
evolution d{F2 n / F2]/ d\h\Q 2 ] provides an important clue to the recombination 
mechanism. 
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